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Abstract 



We derive the convective terms in the damping which determine the struc- 

^ ■ ture of the moving charge-density wave (CDW), and study the effect of a 

current flowing transverse to conducting chains on the CDW dynamics along 

'^ i the chains. In contrast to a recent prediction we find that the effect is or- 

5 ' ders of magnitude smaller, and that contributions from transverse currents of 

O ' 

Q ' electron- and hole-like quasiparticles to the force exerted on the CDW along 

the chains act in the opposite directions. We discuss recent experimental ver- 

■^ . ification of the effect and demonstrate experimentally that geometry effects 

JC , might mimic the transverse current effect. 

cn 

o 

0^ 
0^ 
-(— > 

a 

o 
o 



PACS numbers: 71.45.Lr, 72.15.-v, 72.20.Pa 



Typeset using REVT^ 



Nonequilibrium steady states of driven periodic systems with quenched disorder have 
attracted significant interest. Examples include quasi-lD charge-density waves (CDW)a, 
vortex lattices in superconductoraHl, two-dimensional Wigner crystalsu. All these systems 
are characterized by pinning, so that macroscopic motion is only set up when the applied 
force is larger than a threshold value. 

In the particular case of quasi- ID CDW where the CDW motion is only along the chains, 
many interesting effects were predictedoia such as a dynamical phase transition between a 
disordered and an ordered state and instability due to the proliferation of phase slips. While 
a dynamic transition was observed in the recent experiment on high-quality NbSea crystals^ 
between temporally ordered CDW creep to high-velocity sliding, many of the predicted 
effects have not yet been detected. The predictions were based on the equation of motion 
including a phenomenological convective term in the damping, '-^[dt'-p + (vV)(^], where v 
stands for the velocity of the moving lattice (in our case for the CDW velocity). 

Recently Radzihovsky and TonerB predicted a dramatic effect of transverse current on 
CDW motion along the conducting chains. This effect results in a "current-effect transistor" 
in which the CDW channel is turned on by a transverse quasiparticle current. The effect, 
again, stemmed from a convective term phenomenologically inserted in the equation of CDW 
motion, the velocity v being assumed to be of the order of the normal electron velocity. An 
experimental verification of this effect was reported recently by Markovic et aa. 

The convective term introduced in ref. [^ has a different origin from that of refs. 0-^. In 
the former case such a term reflects the action of the transverse current on the CDW, in 
the latter case it appeared because Galilean invariance was implicitly assumed for damping, 
and the time derivative dt was replaced by the convective derivative dt + vdx- This would 
be valid if damping in the laboratory frame were the same as in the frame moving with the 
CDW. CDW damping is induced by scattering of quasiparticles, thus the scattering would 
be the same in both frames if the scatterers (impurities and phonons) were moving with 
the CDW. Since it is not so, calculation of the damping coefficient from the microscopic 
theory is needed. The explicite calculationli2l yields neither Galilean invariance nor the con- 
vective term. Nevertheless, terms with such a symmetry are allowed, and terms resembling 
the convective one appear for a completely different physical reason, namely, due to local 
perturbations of the quasiparticle density. We discuss this problem shortly below. 

The main purpose of our paper is to elucidate the physical origin and to derive the 
convective term by means of the microscopic transport theoryll2l. We find that this term is 
orders of magnitude smaller than it was proposed in ref. |^; furthermore, it is proportional 
to the difference of transverse electron and hole currents rather than to the total transverse 
current. Thus the net force acting on the CDW is determined by the difference of quasipar- 
ticle contributions from energies above and below the Fermi surface, like in the Hall effect or 
thermopower. In principle, predictions of ref. ^ survive, but the effects become much smaller. 
As for the experimental verification of this effecta, we show that a similar behavior may be 
caused by quite large longitudinal current induced by the transversely applied voltage due 
to the large anisotropy of CDW compounds. 

Though we study here the CDW only, we expect that our conclusions are qualitatively 
valid for other types of moving periodic structures. In these systems the damping is related to 
quasiparticle scattering as well, and the Galilean invariance does not hold, so the longitudinal 
convective term must be proved by a microscopic theory. As for the transverse effect in 



related systems, we expect that it must also be related to the difference between the electron 
and hole current rather than to the total current, as well. 

To derive the equation of motion for the CDW phase we start from the equation for the 
nonequilibrium quasiclassical Green's functions in Keldysh representationHI which we use in 
discrete form with respect to chain indices". The Green's functions 



are formed by the retarded and advanced Green's functions and by g^ which is responsible 
for the electron distribution. Each of g"^ is a matrix with respect to the index identifying the 
sheets of the Fermi surface of the quasi-one-dimensional conductor at +pf and —pr, and to 
the chain number. They satisfy the equation (in units with e = 1, h = 1) 

IVpOxgnm ~r / , <^±i\^nn+i\'^l)9n+iin 9nm+i^m-\-im\J'2)\ 
i 

~r il^yZi J gnm gnmi^m\J'2j ~r ^ng-nm gnm^m U, \^1) 

where matrix product and the convolution with respect to time in the last terms presenting 
the elastic colhsion integral are assumed, /i„ = iazOt + icTyAn — ^nCz, ^n = (t>n — {vF/'2)dx(Pn, 
(pn is the electric potential, A„ and (pn are the CDW amplitude and phase, S„ = ^i^f(Jzgnn<^z — 
j^bWxgnn<^x + crydnnO'y], Vf and v^ are the forward- and back-scattering rates. The term with 
tj_ describes the interchain coupling. A^m = c^ cos ^^^ + i sin ^i^^ '■Pn,m = V^n — V'm- 

Equation for the phase can be found from the self-consistency conditions for the order 
parameterHyO 

Tra.,„„,(M)=.--^^. (2) 

We solve equation (|I|) perturbatively with respect to t±. The function g^^^ we present as 
a sum of regular and anomalous parts, the latter is related to nonequilibrium perturbations 
of the quasiparticle distribution function, g^^ = gnnfm — gnrrf'n + dnm^ ^n = 1 — 2r2^ {e — fj^n), 
Up is the Fermi distribution function, /i„ is the shift of the chemical potential in chain n, 
and e is the energy related to ti — t2 by Fourier transformation. 

In the zeroth approximation we neglect t± and findEl'El g^^^^ = gnC^z + fnio'y, gnn = 

with gn = en/in, fn = ^n/^n, €„ = £-$„ + iu+g^^"^^ /2, U+ = Uf + Uf,, A„ = A„ - 



In the first approximation we calculate Green's functions off-diagonal with respect to 
layer indices. 






the result for retarded (advanced) function can be obtained from (|^) by substituting (t^ 
tn+i) by 1, and all indices A for R {R for A). 

Using (^) we can calculate the transverse current. 

i K f 

jn,n+l{x) = 7T— / t_LiTl{Ann+ign+ln ~ 9nn+l^n+ln) de, 



where k is the inverse Thomas-Fermi screening radius. In the hnear approximation with 
respect to the driving field, E^ = (Vn — Vn+i)/d, related to the difference of electrochemical 
potentials, V^ = $„ — /i„, we get for the simplest case of identical chains and temperatures 

r< A 

jn,n+l{x) = (T±E±, (4) 

ax = a^^- (^cos ^- + ^ j e / smh -, 

where u^ = Uf + z/f,, z/ = z/j + Ub, a±N- is the normal-state conductivity. Under typical 
experimental conditions the phase in the transverse directions varies over large distances 
and one can substitute (pn,n+i for zero. 

If some chains are not in the Peierls state, like in NbSes, an additional contribution to 
the conductivity appears, which is related to the hopping between the chains in the metallic 
state and exhibits no thermal activation. 

The second order corrections in t± to (jnn we calculate by inserting the first order correc- 
tion (^ into equation (|T]) for n = m. After some algebra we get the component of Green's 
function g^ = Tr^^ needed to calculate the equation for the phase by means of eq.(0). In 
addition to standard terms describing the elastic deformation, action of the electric field 



and damping of the CDW which can be found in refs. 0,0, this function contains the term 
which is related to perturbations of the quasiparticle distribution function by the transverse 
quasiparticle fiow. This additional term reads 

5-2 = Yl '^t'iii^n+i - Vn) smipn,n+iG{e), (5) 



G = ^Re 



de^^^ (e+a(e+e„\.)(e^+e^ 



n+iJ 



Integration of expression (|^) over energy gives the force induced by the transverse currents. 
Note that in G{e) the derivative of the Fermi distribution functions is multiplied by an odd 
function of energy. This means that the contributions of electrons (e > 0) and holes (e < 0) 
have different signs. If the chemical potential is not shifted from the midgap position, i. e. 
/i = 0, then the electron and hole contributions compensate each other, and the convective 
term disappears. If /i 7^ then this compensation is not complete. For brevity we assume the 
condition T <^ A which typically holds in CDW conductors practically at all temperatures 
below the fluctuation region near the Peierls transition. Then the equation for the phase 
acquires the form 

m dt"^ dt ^ dx^ 

+2 51 ^li[l + ^(^" ~ Vn+i)] sin ipn,n+i = Fpin + 2VFEn, 
i 

where A = {An/2) J G{e)de, and Fpj„ is a pinning force. In the continuous approximation 
equation the last term in the left-hand side of (0) reads v1_[V']_ip + 2A(E^'V ±)ip], where v'j_ = 
{t±d)'^/2, d is the interchain distance. The term with A is the convective term describing 
the transverse effect. Coefficients 7 and A in the equation depend on details of the energy 




structure of the compound, e. g. on energy dependence of the density of states etc. Noting 
that the main contribution to the integrals stems from the energies close to the gap edge, 
we simplify (^) in two limits, p = (Az/^)^/^/T ^ 1 and p ^ 1. With the accuracy of a factor 
of the order of unity we get 

Both the damping factor 7 and the convective term exhibit a thermally activated behavior 
because they are related to excitations via the Peierls gap. If some chains are not in the 
Peierls state, non-thermally activated contributions appear in the next order approximations 
in t±/A. This happens due to proximity effect which induces non-zero density-of-states 
inside the gap at the chains in the Peierls state and to CDW correlations at the chains in 
the metallic state. Such contributions are of the order of v\;t\l t^ to factor 7, and of the 
order of t5_/^/A^ to coefficient A. Though these contributions are small they may become 
dominant since contributions (^ are thermally activated and may be very small. 

Eq. (|^) shows that the local value of the damping coefficient depends on a local variation 
of the chemical potential, the latter being related to a local perturbation of the quasiparticle 
density bnjn 

7 = 7o(l + a5n/n), (8) 

where a is a factor of order 1 . The last term can be _presented in the form of a convective 
term with the velocity parallel to the CDW motionoEl. We present a "quick and dirty" 
qualitative derivation of such a term. Let us consider the regime of linear screening when the 
shift of the chemical potential b[i from the equilibrium position is small. Then we estimate 
bnjn ~ ^[i-jT in a semiconducting CDW state, and bnjn ~ b^ijEp in the metallic CDW 
state, where Ep is the Fermi energy. b[i can be found from the quasineutrality conditionE3. 
The local perturbation of the charge density consists of the CDW and the quasiparticle 
contribution, p oc [NcVFdx^/2 + (1 — A^c)^/^], where N^ is the fraction of electrons condenced 
in the CDW. Then, expressing 5 11 from the condition p = and inserting it to (|]), we find the 
renormalisation of the damping term. Using the standard expression for the CDW velocity 
V = dt'~p/2kp we present the damping term for the metallic CDW in the form resembling the 
convective derivative 

7 = 7o [dt^ + ^ ^_''^ vd^Vj ■ (9) 

The "convective" term here is proportional to the temperature dependent ratio between 
contribution of the condensed electrons and quasiparticles, which is small in NbSes or m-TaSa 
for the upper CDW region. In the case of semiconducting CDW the additional factor Ep/T 
in the "convective" term appears, and since typically (1 — iVc) -C 1 the renormalisation of the 
damping coefficient is large. Note, however, that in refs. ^|-|^ the convective term contains 
an averaged CDW velocity, while in Eq.@ v is the time and position dependent local CDW 
velocity which can be treated as a constant parameter only in the limit of large CDW velocity. 



Thus electric fields induced by local perturbations of the quasiparticle density renormalize 
the damping similarly to renormalization of the quasiparticle conductivity @ and of the 
CDW stiffnesaly. A detailed discussion of the effects induced by this convective term is 
beyond the scope of the present paper, therefore, we limit the discussion by a comment that 
some predictions of refs. ^|-|^ must be revised. 

Now we return to the transverse convective term. Its sign depends on the sign of the 
charge of the quasiparticles which determine the single-particle conductivity, and its tem- 
perature dependence is determined by an interplay of a thermally activated behavior and 
higher order terms in the interchain interaction. The exact value of the convective term 
depends on the details of the electron spectrum of the material. Using calculated values 
for the factors in equation of motion (|^ one can estimate the magnitude of the transverse 
electric field which may affect the CDW dynamics along the chains. In particular, following 
the calculations of ref. |^ one gets a decrease of the threshold depinning field along the chains 
as a function of the perpendicular electric field E^ at E^ ^ E^, where E^ is the crossover 
field 



£,(£J = S,(0)|e--^/-.. E.^^q^. (10) 

For the case of a semimetallic CDW material like NbSca when non-thermoactivated contri- 
bution to the factor A dominates, the crossover field can be estimated as 

Ee ^ ^ -JEaEt(0), Ea = 10^ ^ lO^V/cm. 

The crossover field Ec is much larger than the longitudinal threshold field Et{0) because 
both Ea and the prefactor before the square root are large. 

The experimental observation of the effect of a transverse current on the longitudinal 
CDW dynamics was reported recently by Markovic et alB in NbSes. It was found that 
application of a transverse electric field of the order of 1 V/cm was enough for a substantial 
reduction of the longitudinal threshold field. The discrepancy with the above estimate is thus 
4-5 orders of magnitude. Now we discuss the problems which may occur in the experimental 
study of the effect in the simple geometry with two side contacts for the application of the 
transverse currentH. Let us consider a sample of width w along y axis with side current 
contacts of length / along the x axis parallel to the chain direction. Assuming a uniform 
current density in the contacts and Ohmic conductivities a\\ ^ cr^ in x and y directions 
(y = at the contact), respectively, we calculate the distribution of the electric potential in 
the sample (see ref. |12|) 

1/2 



■■ Eq / dxi 2_^ In ^ =, (11) 

_//2 "=-°° V (^ - xiyai_ + {y- 2nw) V|| 

where Eq = j±/{2TTa±), j± is the transverse current density through the contact. 

In the case of large conductivity anisotropy, large components of the longitudinal current 
appear near the contacts. According to (pA]) in the region near the contacts, at distances of 
the order of / along the chains and at distances about lJ(y\\/cr± in the transverse direction, 

6 



the component of the electric field parallel to the chains is almost constant, E\\ ^ E{ 
larger distances the electric field decays as 2E\\l/x at distances from x ^ l,y^^/a\\/a_L 
as 2E\\lx/y'^ a.t x,l <^ y 



0. At 
and 



'cr\\/a±. And, finally, it decays rapidly at x > wJa\\/a±. 
Thus, if a transverse voltage V applied to the sample is large enough, E\\ may be larger 
than the threshold depinning field Et- Moreover, the region where E\\ > Ex increases with 
increasing V . Note also that the potential difference applied at different directions from 
the contacts along the chains is equal to V/2: the potential decays along the x axis at 
distances of the order of wJa\\/a±_ from the contact (c/. ([IT])). This voltage may easily 
exceed the threshold voltage needed to drive the CDW along the chains. In the experiments 
of ref. ^ / = 100 /xm, w = 36 /im, <yj^/<y\\ ~ 0.04. Then for the characteristic value of the 
transverse current ~ 200 /iA shown to suppress the longitudinal threshold field, we estimate 
E\\ ~ 0.3 V/cm. That is well enough for depinning the CDW for which Et ~ 0.2 V/cm 



even if one takes into account the phase-slip voltagdl^ of the order of 1 mV required for 
CDW breaking. So we conclude that the effect of the transverse electric field on the CDW 
dynamics can be easily masked by the CDW depinning by the longitudinal electric field 
component which is quite large due to large anisotropy of the conductivity. Moreover, when 
an additional longitudinal current ladd is applied to the central segment of the sample, the 
dependence of the threshold field on ladd niimics that of the transverse current observed 
in ref. ^ (see Fig. p. It also shows the crossover current for the onset of the decrease of 
Eriladd) considered in ref. ^ as an argument against the possibility that the threshold field 
reduction is due to current inhomogeneities around the transverse contacts. The crossover 
behavior may be caused by the phase-slip voltage in the central segment where the currents 
I and ladd have the same direction (see inset in Fig. |I]). 
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FIG. 1. Dependence of the threshold field of segment 2-6 on an additional longitudinal 
current in a geometry shown in inset for two samples of NbSes with contact separations 
h3,k4,k5,k6 » 100 ^m: V - ^ = 115 K, A - T = 128 . o and • - dependence ET{Iy)/ET{0) 
on the transversal current ly (data of ref. 0). ladd is plotted in units It for our data, and ly is 
multiplied by 6 x 10~^ cm^/A for data of ref. 0. 



In conclusion, the calculated convective term is found to be much smaller than it was 



originally suggestedtl. Our analysis of the experimental verification of the effects shows that 
the observed behavior may be caused by geometric effects due to high anisotropy of the 
conductivity. So new experimental studies are needed. 
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